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The inverse- square interacting particle system []I|, |^ and its spin generaliza- 
tions 0, |, H, 13, H, H are important models of many-body systems, due to their 
exact solvability and intimate connection to spin chain systems |TT], |r^, |TB 



2-dimensional Yang-Mills theories |14, 16 1 etc. 



Most of the variants of these systems can be though of as appropriate 'foldings' 
of the basic Calogero model with an augmented number of particles. Versions of 
this idea have appeared in the early literature, and have been been used, e.g., to 
motivate the Sutherland and elliptic (Weierstrass) versions of these systems. 
In this paper, we use this approach in the case of spin-generalized systems to give a 
more intuitive derivation and interpretation to some recently produced models and 
to derive new models. 

We begin with a brief review of known cases. We consider the Calogero model 
in connection with some of its discrete symmetries D. The equations of motion 
remain invariant under the phase space mapping where are phase 

space variables. Then the reduction to the invariant subspace = -D(0) is kine- 
matically preserved, that is, the equations of motion do not move the system out of 
this subspace. Therefore, reducing the initial value data to this subspace trivially 
produces a system as solvable as the original one. The motion will be generated by 
the original hamiltonian on the reduced space. 

The starting point will be the inverse-square scattering particle system 

(1) 

1=1 ^ ^ iytj -^ij 

where Xj, Pi are particle coordinates and momenta in one dimension and 
The considered symmetries are: 

• Translation invariance T: Xj — Xj + a, pi pi 

• Parity P: Xi — > -Xi, pi -pi 

• Permutation symmetry M: Xj — > XM{i), Pi — ^ PM{i) with M any element of the 
permutation group of particles. 

Other symmetries will not be useful for our purposes. 

A direct reduction of the system by any of the above symmetries does not produce 
anything nontrivial or sensible: = T(0) is possible only in the trivial case a = 0, 
while = -P(0) and = M(0) requires (some) of the particle coordinates to coincide, 
which is excluded by the infinite two-body potential. We get useful systems only 
when reducing through appropriate products of the above symmetries. These are: 

D = PM: We reduce by P and a particular permutation: M[i) = N — i + 1 (or 
any other in the same conjugacy class). This is uniquely fixed by the requirement 
that = 1 (so that -D^(0) = M^(0) does not make any two different particle 
coordinates to coincide) and that M{i) = i for at most one i (so that no two or 
more particle coordinates are put to zero). The constraint: 

Xi = -XN-i+l , Pi = -pN-i+1 (2) 
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effectively reduces the original system into two mirror-images. The reduced hamil- 
tonian is 



i=i 2 2 ^ , . {xi XjY 2 ^ , . {xi + XjY 



N' = 



Tl ' ^'=^4 ' -12 
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^' = ^f7 + 2(^M (3) 



where an overall factor of 2 has been discarded and [.], {.} denote integer and 
fractional part, respectively. The second term in the potential is the interaction of 
each particle with the mirror image of each other particle; the third part accounts 
for the interaction of each particle with the mirror image of itself, and with a particle 
fixed at the origin by the constraint (for odd A^). 

We further note that parity symmetry persists in the case where an external 
harmonic oscillator potential is added to the system, promoting it to the confining, 
rather than scattering, Calogero model: 

h-eIpi+It^^+t.^'-^ (4) 

i=l ^ ^ i^j -^ij i=\ ^ 

Reduction by D = PM produces an integrable system similar as above with the 
added harmonic oscillator potential. 

b) D = TM: No finite-rank element of the permutation group will do, since 
repeated application of D would eventually lead to Xi = Xi^ ma. We overcome this 
by starting with N' particles and taking the hmit A^' — > oo. We pick the element: 
M{i) — i-\-N for some finite N which, for infinite N' ^ is infinite-rank. The constraint: 

Xi+N a , pi+N = Pi (5) 

leads to a system consisting of infinitely many copies of a finite system displaced by 
multiples of a. We can parametrize the particle indices by the pair (i,m), i — 1 . . . N 
and m e Z, where the original index is i -|- mN. The constraint now reads 

Xi,m = Xi + ma. , Pi,m = Pi (6) 

The resulting system is infinite copies of an A'"-body system. The reduced hamilto- 
nian is 




In the above summation terms with {i = j,m = n) are omitted since they correspond 
to self-interactions of particles that are excluded from the original Calogero model. 
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The summation over m' above accounts for the infinite periodically repeating copies 
of the system and can be dropped. The infinite m-summation accounts for the 
interaction of each particle with the multiple images of each other particle and can 
be performed explicitly. We eventually get 

1=1 i^j a 

that is, the Sutherland model. In the above we omitted an irrelevant constant term 

2 

equal to coming from terms with i — j, which account for the interaction of 

each particle with its own infinite images. 

c) We can formally extend the T symmetry to complex parameter a. As long 

as there is a subset of coordinates in the reduced phase space that remains real 
and generates all other coordinates through use of we will have a well-defined 
real subsystem. Applying D = TM for infinitely many particles parametrized by a 
double index i,j G Z, for two complex translations a and b, the constraint is 



X. 



^i+N',j + a — Xij+M' + b , Pij — Pi+N',j — Pi,j+M' (9) 



and we end up with a finite system with = M'N' particles periodically repeating 
on the complex plane. Similarly to the Sutherland case, the hamiltonian becomes 
infinitely many copies of 

The above sum has a logarithmic ambiguity that is easily regulated by subtracting 
the constant g'/(ma + nb)^ from each term. We end up with 

N' 1 1 

H = Y.-/.+-^Y.9V{x,,\3M (11) 

i=l i^^j 

that is, the elliptic model with Weierstrass potential. 

d) Di = PM together with D2 = TM: This is a combination of (a) and (b) 
above. We work again with an infinite number of particles. We impose two con- 
straints: 

Xj+N' = + a , Pj+N' = Pj (12) 

where e = 0, 1 (any other choice of e is equivalent to one of these). Parametrizing 
j — i + mN' by the pair i, m, i = 1 . . . N', m e Z, we have 

^N-i+l-e Xi^rn : VN-i+l-e — ~Pi 

^i,m — ''r ma. , Pi,m — Pi (13) 
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so we end up with a finite system of = [(A^' — e)/2] particles. Tlie reduced 
liamiltonian is infinitely many copies of 



2"^ 2 ^ a2 sin^ vr^ 2 ^ sin^ vr^^ 

2 Va^sin^TT^ 2 Va^sin^vrS- ^ ^ 

« a I a 

with 

The coordinates Xi can all be taken in the interval (0, a/2) and the particles interact 
with their infinite mirror- images with respect to mirrors placed at x = and x = a/2 
and with particles fixed at x = (if e = 1) and at x = a/2 (if 2{(A^ — e)/2} = 1). A 
similar construction can be performed with two complex translations, as in (c) plus 
one parity reversal. We obtain a similar model but with elliptic functions appearing 
instead of inverse sine squares. 

The above is fairly standard and exhausts the possibilities for spinless particles. 
Before we proceed to the more novel and interesting case of particles with spin, 
we find it instructive to demonstrate how the above construction reproduces the 
conserved integrals of motion of the reduced system. We will consider case (b), as 
the most generic, case (a) being rather trivial. 



The Lax matrix of the original scattering Calogero model is |T7 

U 

Lij = pAj + (1 - — (16) 

where l"^ = g. Traces of powers of L produce the integrals of motion in involution 
for the model: 

4 = trL\ k = l,...N (17) 

For the system of case (b), we promote the index i into a pair (z, m) and choose 
Xi,m = Xi + ma as in (^). The resulting infinite-dimensional matrix Li^jn can be 
thought of as consisting of infinitely many blocks of size N x N, 171,71 labeling the 
blocks and i,j the elements of each block: 

a 

Lim,jn Pi^mn^ij ~l~ (1 ^mn^ij) | 7 \ l-^^) 

Xij + (m — n)a 

We observe a block 'translational invariance' of the matrix L in the indices m, n, 
which reflects the invariance of the model under a translation by a. Due to this, we 
can trade the pair m, n for a single index m — n 

I'imjn — Lm—n;ij (-^'^) 
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and thus L becomes an infinite collection oi N x N matrices L„ labeled by n. We 
define the Fourier transform L{a): 

L{a)=J2e'^''Ln (20) 

n 

in terms of which L„ is 

^ ^"L((T)e-^"" (21) 



27r JO 

The corresponding integrals of motion are traces of powers of L. Denoting by Tr 
the trace in the infinite-dimensional space labeled by i, m and by tr the trace in the 
A^- dimensional space labeled by i alone, we have: 



Ik — TrL'^ — ^ tr(L„^_„2 ■ ■ ■ Ln^.-ni) 

n\,...n}^ 

{rrii = n-i — rii^i) = ^ ^ tr(Lmi ■ ■ ■ -Z^m^.i-^-mi — m^.i) (22) 

ni mi,...ms;_i 

The sum over rii above produces a trivial infinity. This is due to the summation 
over the infinite copies of the system, just as in the case of the hamiltonian, and 
will be dropped. In terms of the Fourier transformed L[a) the reduced Ik become 
simply 

h = ^ r trL{a)''da (23) 
Zn Jo 

It is now a matter of calculating L{a). From the Fourier transform 

pincr i{n—a)x 

y = — for < a < 27r (24) 

„ + X sin na 

we obtain for Lfcr) 

ini „TT — (y 



'J 



pAj + (1 - — : ^ - i 6i 



(25) 



where the last diagonal term linear in vr — a came from terms with i = j, n ^ in 
Ln,ij- We observe that the matrix inside the square bracket apart from this linear 
part is the standard Lax matrix L of the Sutherland model: 

~ t7r£ 

Lij = PiSij + (1 - Sij) — : ^ (26) 

a 



Substituting (^) in (^) we note that the exponential factors cancel (due to 
Xjjj2 + ■ ■ ■ ^ik-ih = 0) are left with 

Ik = — r%rfL-f^^y^a = Vtr I'^-^iz^ r^Ti-ayda 

2nJo \ a J f^^ s\{k-s)\ 271 a.' Jo ^ ' 



n=0 



(2n+ 1)!(A; - 2n)! V a 



where Ik = iiL^ are the conserved integrals of the Sutherland model. We obtain 
a linear combination of the intergal Ik and lower integrals of the same parity. The 
appearence of the lower integrals originates from the interaction of each particle with 
its own infinite images. We saw an example of such a term in the constant potential 
term that we omitted from the reduced hamiltonian of case (b). In conclusion, we 
have recovered the integrals of the Sutherland model. 

We extend now these considerations to systems of particles with internal classical 
U (n) degrees of freedom. The corresponding starting spin-Calogero system can be 
obtained, for instance, from the model in Q (which can itself be obtained as 
a reduction of a hermitial matrix model |]18[ into nontrivial angular momentum 
sectors) by redistributing the global U{M) degrees of freedom of this model into 
individual particle 'spins'. Equivalently, we can take the infinite-volume classical 



limit of the spin model derived and solved in [15|. The hamiltonian reads 



H = j:lpf + lT.'-^ (28) 

i=l ^ ^ ij^j -^ij 

The Si are a set of independent classical U{n) spins of rank one and length £, that 
is, n X n rank-one hermitian matrices satisfying 

tT{S,y = f (29) 

and with Poisson brackets 

{iSi)ab, iSj)cd} = -iSij [iSi)ad5cb - 5adiSi)cb] (30) 



Such spins can be realized in terms of oscillators [p!5| : 

{Si)ab = A'iAl , a,b = l...n (31) 

where {A°; , A";) are a set of nN independent classical harmonic oscillator canonical 
pairs with Poisson brackets: 

{AlA'^} = i5,,5ak (32) 

and satisfying the constraint 

J2AtA'^ = i foralH (33) 

a 

The above model, in addition to the previous symmetries T, P and M, also 
possesses the symmetry 

• Spin rotations U: Si USiU^^, with U a constant unitary n x n matrix. 
{xi,Pi) remain unchanged. 

Again, reduction by this symmetry alone leads to no interesting system (implying 
either U = 1 or Si = 0). Reduction by PUM or TUM, however, much along the 
lines of the previous PM and TM reductions, produces new and nontrivial results: 
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e) D = PUM with P and M as in (a) before, and U a unitary matrix satisfying 
f/^ = 1 (this is necessary since P and M are of rank two). The constraints are 

Xi = — xtv-i+i , Pi = —pN-i+1 , Si = U S]\f+i_iU ^ (34) 

The reduced hamiltonian acquires the form: 

H = j2-p' + 1^ tr(g,g,) ^ l^ tTjS.US.U-^) 
i=i 2 2 .^^ (xj — Xj)"^ 2 (xj + Xj)"^ 

^ tijS.US.U-') , ^ tr(g,go) 

where A^' = [A^/2], g = 2{A^/2} and So = USoU~^ is an extra spin degree of freedom. 

The form of the hamiltonian for the reduced model and its physical interpretation 
simplifies with an appropriate choice of basis for the spins: by using the [/-invariance 
of the full model, we can perform a unitary rotation V to all spins Si —>■ VSiV~^. 
This transforms the matrix U appearing in (|36| ) into U —>■ V~^UV. With an appro- 
priate choice of V we can always choose U to be diagonal: U = diag{e'"^''). Because 
of the constraint, this means that U will have the form U = diag{l, ...1,-1,... — 1) 
with rii, n2 entries equal to 1, resp. —1. So we see that the original U{n) invariance 
of the model has been broken to U{ni) x U{n2). If ni = n2 there is an additional 
Z2 exchange symmetry. 

As in the spinless case (a), we could have started with a spin-Calogero model in 
an external oscillator potential (which shares the same U and P symmetries), and 
obtain a model as above with the extra confining harmonic potential. 

i) D = TUM with T and M as in (b) before, and U any unitary matrix. The 
constraint on the phase space is 

Xi = Xi+N + a , pi= Pi+N , Si+N = USiU'^ (36) 

The system becomes, again, infinite copies of a-translated and ^/-rotated systems, 
and the reduced hamiltonian is 

In the above we cannot drop terms with i = j any more, since they are now spin- 
dependent rather than constant. Only the term (i = j, m = 0) must be dropped 
from the summation as before. 

Again, the form of the hamiltonian for the reduced model and its physical inter- 
pretation simplifies with an appropriate choice of basis for the spins which makes U 
diagonal: U = diag{e^'^'^). The trace in ( |T7D then becomes 

n 

tiiS.U'^S.U-^) = ^ iS,U{S,),^e-''^'^^'> (38) 

a,b=l 
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where (pat = (pa — (pb- The m-summation appearing in (^) gives 

X! ^ ^"^^2 ^ 5Z {Sdab{Sj)ba X! 7 ^ \2 ^"^^^ 



ni=— oo 



E {S^)ab{S,)ba (40) 



a,fe=l 



with the potential Vab{x) being 



We must distinguish between the cases i ^ j and i = j. For the case ? 7^ j the sum 
can be obtained from the x-derivative of (P^): 

Kfe(x) = \e-'^^^' - ^^<Pab cot — - 7r\(Pab\] (42) 

ysm ^ sl J 

For the case i = j we must omit the term m = from the summation. We obtain 
an x-independent potential: 

V-..>i.(M.)-^)^i-^ (43) 

In the above we omitted a constant (a, 6-independent) term equal to ^ which would 
contribute to the hamiltonian a term proportional to J2i tr(5'i)^. Due to (pOf), this is 
an irrelevant constant. With the above, the reduced hamiltonian eventually becomes 

^11 1 ~ 

^ = E + 2 E E Vabix,,) iS,)US,)ba + 2 E E iS.^ab{S.)ba (44) 

i=l i^j a,b i a,b 

This is a model of particles with U (n) spins interacting through U (n) non-invariant 
couplings, due to the presence of the matrix Vab- The original global U (n) invariance 
is, now, broken to the diagonal f/(l)" part and only the diagonal components Saa 
of the total spin 

Sab = J2(S^)ab (45) 

i 

are conserved. The standard [/(n)-invariant spin-Sutherland model is recovered 
upon choosing Vab ~ Sab, in which case the sums over a, b above become a normal 
trace. This is achieved by choosing const ant, that is, U = e**^. 

The above model is, in fact, the same as the classical model introduced by Blom 
and Langmann [rU|, and this author [^, in the particle-spin form in which it was 



recast in 20 



H = lT.P' + lT.{T.ya^i^^MU{S^^^^ (46) 

t i^j \ tj 2 J 

+ ^YU2^^b{Si)aa{Si)bb + ^J2^'^b{qaqb- SaaSbb) (47) 
i ab ab 



To fully see the equivalence, we must observe the following: 

1. In the present construction we expressed the hamiltonian in terms of U{n) 
spins Si. In it was, instead, expressed in terms of traceless SU{n) spins Si. By 
(pip and (|33|) we have trS'j = i, so the relation between the two is 



S^ = S,-- (4J 



n 



2. The expression (^) derived in |^ was fully quantum mechanical. It can 
be seen that the term i{i + n) in (|47|) classically becomes i"^ {n was a quantum 
correction similar to the shift of the classical angular momentum to J{J + 1)). 

3. For the rank-one matrices Si we have the relation 

{Si)ab{Si)ba = {Si)aa{Si)bb (49) 

4. In a set of dynamically conserved charges were introduced that can be 
chosen to have any value as long as they sum to zero. 

5. In 1^ the particles were taken to move on the unit circle, that is, a = 27r. 



Doing the above substitutions in (HI) we see that it becomes practically identical 



to (|47|) . The two expressions differ by constant terms depending on the charges 
and the diagonal elements of the total spin Saa- Since both of these quantities are 
constants of the motion, the two models are trivially related. 

g) We can, similarly to (c), extend the above construction to two complex trans- 
lations and corresponding spin rotations. Parametrizing again the infinite number 
of particles with a doublet of indices j G Z the constraints are 

= Xi+N,j + ^ ; Pi.j = Pi+N,j 5 Sij^Nj = U SijU 

= Xi,j+M + b , Pij = Pijj^M , Sij+M = ^SijV ^ (50) 

and we end up as before with a finite system with N' = MN particles periodically 
repeating on the complex plane. Since the two space translations and the corre- 
sponding particle permutations commute, for consistency the two spin rotations 
must also commute: 

UV = VU (51) 
The corresponding reduced hamiltonian becomes infinite copies of 

Just as in the previous case, we can choose a basis for the spins that diagonalizes 



both U and V to U = diag{(f)i), V = diag{6i). The m, n-sums that appear in (|52D 
become 

n oo ^-im<j)^f,-ineab 
E i.Si)ab{,Sj)ba E 7 \ \ rvl (53) 
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where the term m — n — is omitted Hi — j. We obtain again a potential Vab{xij), 
for i ^ ji given by the sum 

and a spin self-couphng Vab for i — j, given by 

g-im^ai,— • 

(m,n)7^(0,0) 



^«f= E 7 —-13^ (55) 



Note that now, due to the presence of the phase factors, these sums are convergent 
and have no regularization ambiguity. The only ambiguous terms, defined modulo 
an additive constant, are the ones with a — b. We will comment on the impact of 
such a regularization ambiguity in the sequel. 

The potential Vab{x) is a modular function on the complex torus (a, b) with 
quasiperiodicity 

Vab{x + h) = e'''^>'Vab{x) (56) 
It has a double pole at x = 0, with principal part 

K(,(2:) = ^ + 0(x°) (57) 

and no other poles in each cell. These properties uniquely define Vab and allow for 
an expression in terms of theta-functions. We put 

•.^ Oi (^x - q^)) 01 (^{x - q2)) 

where qi^2 are the as yet unknown zeros of Vab{x) and the theta-functions appearing 
above have complex period T — b/a. This has the right quasiperiodicity under 
X ^ X + a.. In order to also have the right quasiperiodicity under x — > a; + b, gi^2 
must satisfy 

?1 + g2 = 7i-(a^a& - b^ab) = Qab (59) 
ZTT 

and to have the right behavior around x = we must further have 



7r 



3? 



(60) 



^1 {iqi) 9, {iq^) ^ 
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The equations and ( pTD above determine qi and while (^) then determines 
A. It may be possible to express qi, q2 in a more explicit form, or to recast (|58D 
in a form more symmetric in a,b, by using theta-function identities. Finally, the 
self-coupling Vat can be extracted from Vab{ 

Vab = lim (Vabix) - ^) (62) 

To sum up, we obtain a U{n) non-invariant spin-generalization of the elliptic 
model given by a hamiltonian of the form (|^) but with the potentials appearing 
now being given by ( |58| , |62D . The U{n) invariance of the original model is, again, 
broken down to the diagonal abelian sungroup [/(l)"" due to the dependence of 
the potential on a, b. The f/(n)-invariant spin-Weierstrass model is regained for 
4>ab = ^ab = 0, that is, trivial matrices U and V. 

We point out that for 9ab = (pab = 0, that is, Qab = 0, the equations for gi 2 
(^,^) are satisfied for any qi = —q2 leading to an apparent arbitrariness. As can 
be seen, however, by applying the addition formula 

9,{x + q)9i{x - q)9,{0f = 9r{xf9,{qf - 9,{xf9,{qf (63) 

this simply amounts to an arbitrary additive constant to the expression for Vab{x) = 
V{x). This corresponds to the need for regularization for this expression in the 
absence of phases, as explained before. (In the case of the Weierstrass function this 
is fixed by further requiring that the 0(x°) part of the function at a; = vanish, 
which picks q = 7tT/2 and makes 94{q) above vanish.) We also point out that we can 
pick any of these values for qi = —q2 at the limit Q = by appropriately choosing 
the ratio (pab/^ab as they both go to zero. 

The ambiguity of the terms with a = b can be fixed in the same way: we can 
choose phases 0aa 7^ 0, 9aa 7^ 0, evaluate the expressions, and then let (paa, (^aa 0. 
This will lead to arbitrary additive constants Ca, depending on the ratio (t)aa/9aa as 
we take them to zero. The same constants, however, will appear in both Vaa{x) and 
Vaa- Their net contribution to the hamiltonian will be 

^H=WT.Ca {S.)aa{S,)aa + ^ E E {S.)aa{S.)aa = ^ E {Saaf (64) 

Since the diagonal components of the total spin S are still constants of the motion, 
due to the residual [/(l)"" invariance, this amounts to the addition of an overall 
constant, and thus leads to systems that are trivially related. 

The same discussion applies if the angles 0^ and 9^ coincide for two or more 
values of a belonging to a subspace of indices J, in which case 0afe = 9ab = ^ for 
a,b & I. This will result to a constant additive matrix Cab in the potential for this 
subspace of indices /, leading to an extra contribution to the hamiltonian 

Aff= E CabSabSba (65) 
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Since (pat = Oat = ^ for a,b G /, however, the corresponding subgroup of U {n) 
remains unbroken, and thus the corresponding components of the total spin Sab 
appearing above are constants of the motion. Once again, the arbitrary terms are 
constant and we essentially obtain a unique system. 

Overall, this is a generalization of the spin-Weierstrass model to one involving 2n 
phases that break the U{n) invariance and promote the potential to a modular func- 
tion. The potential lives on a complex torus in the coordinates, where translations 
around each nontrivial cycle are accompanied by spin transformations. The model 
obtained in (f) can be though of as the limit of the present model with b — > cxd, in 
which case 6a become irrelevant. The properties of this modular potential and of 
the corresponding new integrable model deserve further study. 

We conclude by commenting on the quantum mechanical extensions of the new 
models introduced here. All known classical models of Calogero type are also in- 
tegrable quantum mechanically. We expect, therefore, that the models introduced 
here will also have integrable quantum mechanical counterparts. This is indeed the 
case at least with the model (f), which was fully solved quantum mechanically in 
20]. The quantum version of these models, however, does not seem to be directly 



accessible by the method used here. The reason is that, in general, the constraints 
implied by the restrictions are second class. Therefore, we cannot carry over the 
solution of the quantum model in the unrestricted space and apply the constraints as 
operator relations on the Hilbert space. The proper quantum mechanical treatment 
of these new models remains, therefore, an interesting open issue. 
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